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Recently, a general formalism was presented for gyrotropic, ferroelectric, and multipolar order in 
spin-orbit coupled metals induced by spin-spin interactions. Here, I point out that the resulting order 
parameters are equivalent to expectation values of operators that determine natural circular dichroic 
signals in optical and x-ray absorption. Some general properties of these operator equivalents and 
the resulting dichroisms are mentioned, and I list several material examples in this connection, 
including Weyl semimetals. The particular case of the tensor order in the pyrochlore superconductor 
Cd 2 Re 2 C >7 is treated in more detail, including calculations of the x-ray absorption and circular 
dichroism at the O K edge. 
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Reciprocal (natural) and non-reciprocal dichroism, and 
analogous effects in resonant x-ray scattering, are well 
known probes of materials which break inversion or time- 
reversal symmetry, respectively. There has been a resur¬ 
gence of interest in these effects due to their prominent 
observation in multiferroics and chiral magnets. It has 
been recently realized that similar manifestations occur 
in novel topological materials. For instance, circular 
dichroism due to the chiral anomaly in Weyl semimet¬ 
als has been proposed by Hosur and Li [1] and observed 
by Kerr rotation in Cd 3 As 2 [2]. This is effected by apply¬ 
ing parallel electric and magnetic fields in order to realize 
the chiral E ■ B term, which has been earlier studied in 
the magnetoelectric material Cr 2 0 3 [3]. 

Much of the interesting phenomena associated with 
topological materials can be traced to momentum-spin 
locking due to spin-orbit coupling. Fu has pointed out 
that spin-spin interactions in the presence of spin-orbit 
coupling not only lead to magnetism, but other types of 
order due to this same momentum-spin locking [3]. This 
gives rise to the free-energy 

Ess = ^Fos(A:)•s(fc , ) + ^ 1 l(fc•s(fc))(fc , •s(fc , )) (1) 

k,k' 

+F 2 (k x s(k))(k' x s(k')) + F 3 Qij(k)Qij(k') 
where the quadrupolar term is 

Qijik) = \{.hsj{k) + kjSi{k)) - *0 • s(k))Sij ( 2 ) 

with k the momentum and s the spin. For each term, 
one can specify a particular order - ferromagnetism for 
Fq, isotropic gyrotropic order for F\, ferroelectricity for 
F 2 , and pseudo-deviator order for F 3 - the last has been 
suggested for the pyrochlore superconductor Cd 2 Re 2 07 
[5] . Besides Fo, these all break inversion symmetry, not¬ 
ing that a pseudo-deviator is a rank-2 tensor that is odd 
under inversion. Moreover, as also discussed in Ref. Tj 
each term leads to a unique splitting of the Fermi surface 
(for some illustrations of this splitting, see Ref. [B]) - 


Here, I wish to point out that the resulting operators 
listed above [k ■ s(k ), etc.) are equivalent to operators 
that have been suggested to determine optical and x-ray 
dichroism mm if a mapping is made between k and t , 
where t is the polar toroidal moment mm- k and t have 
the same transformation properties under inversion and 
time reversal (both being odd) [IT], and thus for the ar¬ 
gumentation made in Refs. [7] and [5] they are equivalent. 
This can also be appreciated from the analogy between 
the spin-orbit interaction, (k x VK) • s, and the toroidal 
contribution to the free energy, — t ■ (P x M ), where V 
is the crystal potential, P the polarization, and M the 
magnetization. Polar toroidal ordering has been seen in 
a number of multiferroics mm, and is also the basis 
for a novel theory of the pseudogap phase of cuprates by 
Varma and collaborators ns. following earlier theoret¬ 
ical suggestions of Gorbatsevich et al. QS- In general, 
there are several operator equivalents that can describe 
natural dichroism as listed in Table 1, including a re¬ 
cently advocated form involving the Berry curvature by 
Zhong, Orenstein and Moore [IB] , 

To understand the consequences of this, I discuss each 
order parameter in turn, denoting each as where 
i is the index of F, in Eq. 1 (that is, X 0 = ^ fe s(fc), 
X\ = k ■ s(k), etc.). Xq is ferromagnetism, and can 
be determined by XMCD (x-ray magnetic circular dichro¬ 
ism) [ 21 ]. The rest, though, are associated with natural 
(as opposed to non-reciprocal) optical activity, and cor¬ 
respond to HHE (piezo-electric) terms in the free energy, 
where H and E are the magnetic and electric fields [7]. 
Xi is a pseudo-scalar (also known as an axial toroidal 
monopole | 19 |). It can be seen in natural optical dichro¬ 
ism, due to an interference between electric dipole (El) 
and magnetic dipole (Ml) terms, and the condition for 
its observation for various point groups has been tabu¬ 
lated by Jerphagnon and Chemla 1201 . with a few relevant 
examples shown in Table 2 . Recently, natural dichroism 
has also been seen in connection with spin wave excita¬ 
tions in a field induced chiral magnet, where it is due 
to dynamic terms in the diagonal part of the magneto- 
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TABLE I: Various operator equivalents for natural circular 
dichroism. These are formed from the vectors listed in the 
first two columns, with the transformation properties of these 
vectors under inversion and time-reversal (respectively) listed 
in parenthesis. The scalar product of the two is responsible 
for pseudo-scalar optical activity, and can be shown to vanish 
in the single particle approximation because of the orthog¬ 
onality of the two vectors. The cross product of the two is 
responsible for polar vector optical activity. The remaining 
five components involving the two vectors is a traceless sym¬ 
metric matrix, responsible for pseudo-deviator activity that 
dominates in the x-ray regime. Here, m is the magnetization, 
t the polar toroidal moment, d the electric dipole moment, 
g the axial toroidal moment, v the velocity, Sis the Berry 
curvature, and k the momentum. 


Vector 1 

Vector 2 

Reference 

m (+,-) 

t (t) 

[Hill] 

d (->+) 

g (+>+) 

PUCES] 

v (-,-) 

9s (+r) 

m 

k (t) 

m (+,-) 

this work 


TABLE II: Examples of several non-centrosymmetric mate¬ 
rials and the types of natural optical activity that can occur 
|20j . where PS stands for pseudo-scalar, PV for polar vector, 
and PD for pseudo-deviator (with 1 denoting allowed). Just 
because it is allowed does not necessarily mean dichroism can 
occur. For instance, for Cd 2 Re 207 (I4m2 space group), TaAs, 
and Bi2212, XNCD is zero for x-rays propagating along the c 
axis, due to the presence of glide and mirror planes. 



Space Group 

PS 

PV 

PD 

HgTe 

F43m 

0 

0 

0 

BiTel 

P3ml 

0 

1 

0 

TaAs 

I4imd 

0 

1 

0 

SrSi 2 

P4 3 32 

1 

0 

0 

Cd^R^Oy 

I4i22 

1 

0 

1 


I4m2 

0 

0 

1 

LiOsOa 

R3c 

0 

1 

0 

CdS, ZnO 

P63H1C 

0 

1 

0 

Te 

P3i21 

1 

0 

1 

Bi 2 Sr 2 CaCu 2 08+1 

Bb2b 

0 

1 

1 


electric susceptibility tensor [22] . In general, one expects 
natural dichroism to be present for chiral magnets since 
S x S = iS, leading to a k ■ S term in scattering [23]. In 
the x-ray regime, it is difficult to observe since the Ml 
matrix elements for excitations out of core orbitals are 
typically very small (though they can be significant for 
shallow core levels). 

Moreover, Marri et al. [24] have pointed out that the 
pseudo-scalar vanishes in the single-particle approxima¬ 
tion. This can be appreciated from Table 1, where m 
and t are orthogonal vectors (t being the cross product 
of r and m), meaning that they must refer to different 


electrons to obtain a non-zero result. The same logic [T5] 
also applies to d and g ( g being the cross product of r 
and d) . The same result has recently been pointed out in 
the Berry curvature approach by Zhong et al. Hi , since 
again, v and 9 b are orthogonal vectors (v being the cross 
product of E and fis). The author has verified this using 
the FDMNES program (described below) for tellurium 
(Table 2), where the El-Ml (electric dipole - magnetic 
dipole) interference terms are found to be traceless for 
the various edges studied (Ly L 2 j 3 , M 4j5 , N 4j5 ), noting 
that the E1-E2 (electric dipole - electric quadrupole) in¬ 
terference terms are traceless by definition. Obviously, 
the single-particle approximation cannot be valid in gen¬ 
eral, since natural optical activity can be seen in liquids 
where the pseudo-deviator averages out [25]. 

As commented in Ref. (U X 2 is equivalent to ferroelec¬ 
tric order. The dichroism associated with this is the vec¬ 
tor part of the natural optical activity, commonly known 
as Voigt-Fedorov dichroism. It is due to a longitudinal 
component of the electric polarization, a specific prop¬ 
erty of pyroelectric materials [20]. Typically, the small 
longitudinal electric field is difficult to observe, but one 
can easily study a related effect where the polarization 
of the reflected beam rotates if the propagation vector 
of the incoming beam is not along the surface normal 
[201 . This effect has been seen optically in CdS (Table 
2) [27]. Recently, it has also been seen by resonant x-ray 
scattering in ZnO as an x-ray circular intensity differen¬ 
tial .28]. Since the latter measurements were done at the 
Zn K edge, the resulting operator actually involves the 
orbital and orbital toroidal moments rather than their 
spin counterparts [7]. A related non-reciprocal effect for 
generating a longitudinal electric field in topological Weyl 
semimetals that violate time-reversal symmetry has been 
pointed out by Kargarian et al. [29] ■ As in CdS, a rota¬ 
tion in the polarization in the reflected beam should also 
be seen in this case if the incidence direction is not along 
the surface normal [2S] , 

The non-centrosymmetric point groups associated with 
this order may or may not have a pseudo-scalar compo¬ 
nent as well [20] (Table 2). Note that for toroidic order, 
the application of an electric field perpendicular to the 
toroidal moment induces a magnetic field perpendicular 
to both, and vice versa [30]. A similar ‘inverse Rashba 
effect’ has been discussed in Ref. 22for Rashba spin-orbit 
coupled metals (with k x S along z), where current ap¬ 
plied along y leads to a magnetic field along x (and vice 
versa). For Dresselhaus coupling, one finds a longitudi¬ 
nal response instead [5]. Optical activity can be induced 
as well by currents, the current operator j transforming 
the same way under parity and time-reversal as k and 
t. This has been demonstrated in tellurium, where a 
current induced (non-reciprocal) optical rotation is seen 
in addition to the zero current (natural) rotation due to 
the chiral space group EH- Ref. H] also mentions the re¬ 
cently claimed ferroelectric metal LiOsOa [32]. At the 
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structural transition, the space group converts from the 
centrosymmetric R3c to the non-centrosymmetric R3c. 
According to Ref. [20] the latter space group should give 
rise to vector optical activity, but no pseudo-scalar or 
pseudo-deviator, and in that sense is similar to CdS and 
ZnO, as well as the Weyl semimetal TaAs [33] (Table 2). 

Most interesting is A 3 , which is known in the dichroic 
literature as a pseudo-deviator [ 20 ], and again may or 
may not be observed depending on the space group. This 
operator is largely responsible for x-ray natural circular 
dichroism (XNCD) |7], and is due to electric dipole (El) 
- electric quadrupole (E2) interference resulting from in¬ 
version breaking (circular dichroism being the difference 
in absorption between left and right circularly polarized 
light). At the K edge, it again involves orbital as opposed 
to spin operators. There are five possible terms (since Q 
is a traceless symmetric matrix), and depending on the 
crystal symmetry, they can form one-dimensional, two- 
dimensional or three-dimensional representations. As 
discussed in Ref. |3] a double-dimensional representa¬ 
tion has been proposed for the tetrahedral pyrochlore 
Cd 2 Re 2 C >7 [33] with the associated Goldstone mode seen 
by Raman scattering [35]. Though no optical activity 
has been reported, second harmonic generation has been 
seen [5], 

When considering such tensor order in Cd 2 Re 2 C> 7 , note 
that the double-dimensional representation has as a ba¬ 
sis the two space groups, Mi22 and I4m2 [33] . Below 
the structural ordering temperature, this degeneracy is 
broken, and the latter crystal structure has been claimed 
from x-ray measurements [38] . On the other hand, it has 
been proposed that the former is realized below a second 
lower structural transition temperature m- This can 
be determined by XNCD. Assuming x-rays propagating 
along the c axis of the tetragonally distorted crystal (ad¬ 
vantageous as well, since it eliminates any contamination 
from linear dichroism), the XNCD signal is zero for the 
latter crystal structure due to the presence of a mirror 
plane. But it can exist for the former. 

This is demonstrated by explicit calculations at the O 
K edge with the multiple scattering Greens function code 
FDMNES [38] including spin-orbit interactions [39]. The 
simulations were done using local density (LDA) atomic 
potentials (Hedin-Lundqvist exchange-correlation func¬ 
tion) in a muffin tin approximation that considers multi¬ 
ple scattering of the photoelectron around the absorbing 
site [40]. The cluster radius is limited by the photoelec¬ 
tron lifetime [IT] . Results shown here are for a cluster 
radius of 6 A, though some calculations were performed 
for a radius of 7 A with similar results. Atom positions 
were taken from Huang et al. [32] for the space groups 
I4m2 (which has no dichroism) and I4i22. A Hubbard U 
of 7 eV [42] was used on the Cd sites since the 3d bands 
sit too high in energy in the LDA. Note there are two 
general types of oxygen atoms, six Oi and one O 2 per 
formula unit, Oi having four (three) distinct crystallo- 
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FIG. 1: (Color online) Calculated X-ray absorption (a) and 
circular dichroism (b) near the O K edge for Cd 2 Reo 07 in 
the I4i22 space group for k vectors along three directions - 
(201), (001) and (100) in tetragonal notation. Note that (201) 
corresponds to (111) in cubic notation. For the dominant Oi 
sites, the average cluster Fermi energy is at -6.9 eV. The same 
is shown in (c) and (d), where now the Re sites are assumed 
to be in a d 3 p 2 configuration instead of the d 5 configuration 
used in (a) and (b). In this case for the dominant Oi sites, 
the average cluster Fermi energy is at -8.6 eV. 


graphic sites in I4m2 (I4i22), and O 2 two (one) in I4m2 
(I4i22). 

Results are shown in Fig. 1 for various propagation vec¬ 
tors of the x-rays. The x-ray absorption itself (Fig. la) 
is almost insensitive to the propagation direction since 
the crystal structure is nearly cubic. The dichroism, of 
course, is very sensitive to this (Fig. lb). It is essentially 
zero for a propagation vector of ( 201 ) in tetragonal nota¬ 
tion, corresponding to (111) in cubic notation. This is of 
note, since all experiments referred to here on Cd 2 Re 207 
have been conducted on samples exposing a ( 111 ) face. 
A sizable dichroism, though, is found for (100) and (001) 
vectors (right panel). The absorption at the O K edge 
for Cd 2 Re 2 C >7 has been reported by two groups [321 133] , 
Three prominent features are found in the data, some¬ 
what similar to the calculated features in Fig. la in the 
range of -10 eV to 10 eV, though the lowest energy fea¬ 
ture is not split into two as in the calculation, and the 
middle feature is more peaky in the calculation (the ex¬ 
periment having a broad shoulder extending to higher 
energies). The correspondence to the experimental data 
is improved by putting the Re sites into a 5 d 3 configu¬ 
ration, promoting the other two 5 d, electrons into the 6 p 
states (Fig. lc). The resulting dichroism signal (Fig. Id) 
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is similar to the previous case. The correspondence could 
be further improved by increasing the core hole and pho¬ 
toelectron lifetime broadening [41]. Regardless, an ap¬ 
preciable dichroism is predicted, and if looked for, could 
be used to tell which space group is actually realized be¬ 
low the upper and lower structural phase transitions. 

It may seem puzzling that natural optical activity is 
being described by products of operators that each break 
time reversal. As commented on by Di Matteo et al. UHL 
although XNCD is often claimed to be due to the rank- 
2 tensor product of the orbital angular momentum and 
orbital toroidal moment (row 1 of Table 1) [312], a more 
natural explanation is due to a quadrupolar arrangement 
of electric dipoles, known as an axial toroidal quadrupole 
(row two of Table 1). Such is easily demonstrated to be 
present in V 2 O 3 , for instance, due to movements of the 
oxygen atoms off their high symmetry locations to form 
counter propagating dipole moments on the hexagonal 
planes above and below the vanadium sites [18] . 

Finally, more detailed information can often be ob¬ 
tained from x-ray scattering, since one can obtain infor¬ 
mation as a function of in-coming and out-going polariza¬ 
tion, as well as azimuthal angle. An additional advantage 
is that local inversion-breaking effects can also be seen in 
centrosymmetric space groups due to the differing phase 
factor on each atom that arises from Bragg vectors not 
at the origin Hl| . 

In summary, natural circular dichroism is a sensi¬ 
tive probe that can be used to identify the inversion- 
symmetry breaking order parameters suggested in Ref. 3] 
In addition, if time-reversal is broken, other novel types 
of dichroism can occur as well, including magneto-chiral 
and non-reciprocal linear dichroism [7], along with their 
x-ray scattering analogues m- In general, dichroism and 
related phenomena should be powerful probes of topolog¬ 
ical materials, particularly Weyl semimetals that break 
either inversion or time-reversal symmetry. 
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